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Abstract 

We perform Monte Carlo simulations of a three-dimensional spin system with a Hamil- 
tonian which contains only four-spin interaction term. This system describes random 
surfaces with extrinsic curvature - gonihedric action. We study the anisotropic model 
when the coupling constants Ps for the space-like plaquettes and Pt for the transverse- 
like plaquettes are different. In the two limits jSs = and /St = the system has been 
solved exactly and the main interest is to see what happens when we move away from 
these points towards the isotropic point, where we recover the original model. We find 
that the phase transition is of first order for /3t = /^s ~ 0.25, while away from this point it 
becomes weaker and eventually turns to a crossover. The conclusion which can be drown 
from this result is that the exact solution at the point /J^ = in terms of 2d-Ising model 
should be considered as a good zero order approximation in the description of the system 
also at the isotropic point (3s = (3t and clearly confirms the earlier findings that at the 
isotropic point the original model shows a first order phase transition. 



1 Introduction 



In this article we shall consider a model of two-dimensional random surfaces embedded into 
a Euclidean lattice Z^, where a closed surface is associated with a collection of plaquettes. 
The surfaces may have self-intersections in the form of four plaquettes intersecting on a 
link. Various models of random surfaces built out of plaquettes have been considered in 
the literature |^. The gas of random surfaces defined in corresponds to the partition 
function with Boltzmann weights proportional to the total number of plaquettes. In this 
article we shall consider the so-called gonihedric model with extrinsic curvature action 
0, 0. The gonihedric model of random surfaces corresponds to a statistical system with 
weights proportional to the total number of non-flat edges n2 of the surface The 
weights associated with self-intersections are proportional to kn^ where n^^ is the number 
of edges with four intersecting plaquettes, and k is the self-intersection coupling constant 
[0, Q . The partition function is a sum over two-dimensional surfaces of the type described 
above, embedded in a three-dimensional lattice: 

z{(3)= y: ^-''^""^ (1) 

{surfaces M} 

where e(M) = + 4/cn4 is the energy of the surface M. 

In three dimensions the equivalent spin Hamiltonian is equal to 



Hgonihedric " X! 0-f^O-f+^ + 2 ^ ^^^r+a+/3 7^ X! ^f'^'f'+"^r+(3+/3'^r+/3' C^) 



k ^ 1-k 
and it is an alternative model to the 3D Ising system 



' r'^r+a • 

r,d 



The degeneracy of the vacuum state depends on self- intersection coupling constant k |^. 
If 7^ 0, the degeneracy of the vacuum state is equal to 3 ■ 2^ for the lattice of size N^, 
while it equals 2^^ when k = 0. The last case is a sort of supersymmetric point in the 
space of gonihedric Hamiltonians 

^gonihedric ~ Cr^Cr^+QCr^_|_5+^Cr^_|_^- (3) 

This enhanced symmetry allows the construction of the dual Hamiltonian which has the 
form 

Htai = - E [^"(0 ■ + x) + R'{i) ■ R'{i + ri) + R'iO ■ R'i^ + 0] , (4) 

where V and <j are unit vectors in the orthogonal directions of the dual lattice and R^, 
W and are one-dimensional irreducible representations of the group Z2 x Z2. 

To study statistical and scaling properties of the system one can directly simulate 
surfaces by gluing together plaquettes with the corresponding weight exp{—j3{n2 + 4fcn4)) 
or (much easier) to study the equivalent spin system (Q). The first Monte Carlo simula- 
tions P, 1^, |, 3 demonstrate (see Figure 0) that the gonihedric system with intersection 
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Figure 1: Phase diagram for different values of intersection coupling constant k. 



coupling constant greater than kc ~ 0.5 (including A; = 1), undergoes a second order 
phase transition at f3c ~ 0.44 and that the critical indices are different from those of the 
3D Ising model. Thus they are in different classes of universality. On the contrary, the 
system shows a first order phase transition for k < kc, including the "supersymmetric" 
point k = 0. 

Essential progress in our understanding of the physical behavior of the system has been 
achieved by means of the transfer matrix approach p, |TD[. The corresponding transfer 



matrix can be constructed for all values of the intersection coupling constant k |]rT| and 
describes the propagation of the closed loops in the time direction. In this article we shall 
consider only the k = case. The corresponding transfer matrix has the form [B] 



K{Qu Q2) = exp{-/5 [k{Q{) + 21{Q, A Q^) + k{Q2)] }, 



(5) 



where Qi and Q2 are closed polygon-loops on a two-dimensional lattice, k{Q) is the 
curvature and 1{Q) is the length of the polygon-loop Q This transfer matrix describes 
the propagation of the initial loop Qi to the final loop Q2 and corresponds exactly to the 
Hamiltonian (^,(^). Thus in order to study the critical behavior of the = system 
(|D,(0), one should find the spectrum of the transfer matrix (^). Generally speaking 
all three dimensional problems of statistical mechanics are extremely complicated and in 
our case the exact solution still remains out of reach. However the breakthrough comes 
from the exact solution of a closed system when the transfer matrix depends only on the 
symmetric difference of initial and final loops Qi A Q2 ||10|, |lT| 



A'(gi,Q2)=exp{-2/3/(QiAQ2) }• 



(6) 



The spectrum of the last system has been evaluated analytically in terms of correlation 
functions of the 2D Ising model. This is a nontrivial example of exactly solvable system in 
three dimensions. In particular the largest eigenvalue Aq is exactly equal to the partition 
function of the 2d Ising model and the free energy of our 3d system is therefore equal to the 
free energy of the 2d Ising model. This result nicely explains why the critical temperature 
of the three-dimensional gonihedric system is so close to the critical temperature of the 
two-dimensional Ising model 2(3c = ^ln{l + \/2) ^ 0.44. The next to the largest eigenvalue 

^We shall use the word "loop" for the "polygon-loop" . 
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Ai of the 3d system (|]) coincides with internal energy u{P) of the two-dimensional Ising 
model 

^ = i^r^r+i) = -n{f3) (7) 

and the correlation length defined through the ratio of eigenvalues Ai/Aq, is equal 
to: ^ ^ 

^^^^ " -/n(Ai/Ao) " 

If the internal energy u{j3) approaches 1, the correlation length tends to infinity and 
signals a second order phase transition in the system. But the internal energy of the 2d- 
Ising model drastically increases at the critical point 2/3^ = ^ln{l + a/2) without reaching 
the value 1 ! Thus, in three dimensions, we have the extraordinary situation that the 
specific heat has the logarithmic singularity of the 2d-Ising model, but the correlation 
length remains finite. The conclusion is that the system undergoes a weak first order 
phase transition rather than a second order phase transition. The Hamiltonian which 
corresponds to the transfer matrix has been found in |[ll[] and is equal to 



Hq,aq2 = J2 (8) 

Ex ,Ey 

where a summation is only over vertical plaquettes, the interactions take place only on 
the vertical planes and Ey and the "horizontal" interactions have been switched off. 

The above consideration poses the following interesting question: let us consider the 
system 

Hanisot = PsYl ^^^^ + Pt crcrcrcT, (9) 

Ez Ex,Ey 

which has anisotropic coupling constants for vertical and horizontal plaquettes. We have 
seen that when [3s = ^ the system reduces to the system ® ,(|]) and as we just explained 
it has been solved in [ffO, O] . When /3t = it factors into identical two dimensional plane 



systems solved in and it is always in the disordered phase. Finally when Ps = Pt 
we arrive at our original k = system (|l),(i),(iD- Thus we know the behavior of the 
system at both /^^ = and /3t = 0, but we still don't know the analytical solution at 
the isotropic point Ps = Pt- The understanding of the phase structure of the anisotropic 
system on {i3s,Pt) plane by means of Monte Carlo simulations can drastically clarify 
the situation. Indeed, the important question to which we would like to find an answer is 
whether or not there is any dramatic changes in the behavior of the system when we move 
out of solvable point /^^ = to the isotropic point Ps = Pt, where a first order phase 
transition has been observed in |0, P). Our Monte Carlo simulations show that there are 
no changes in the behavior of the system as it moves from Ps = to the isotropic point 
Ps = Pt- The conclusion which can be drawn from this result is that the exact solution 
at the point Ps = in terms of 2d-Ising model should be considered as a good zero order 
approximation in the description of the system also at the isotropic point Ps = Pt if one 
consider a perturbation with the coupling constant Ps/ Pt- This should be checked by 
further analytical consideration. 
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2 The lattice model 

Thus the lattice action @ may be written in the form 

Pxyif) = 1 — a{r)a{f + x)a{f+ x + y)a{f+ y), 
Pxzif) = 1 — a(r)cr(r + x)a{f+ x + z)a{f+ i), 

Pyzif) = 1 — a{r)a{r + y)a{f + y + z)a{r + z). 

Our goal is to find the phase diagram in the extended (3s ~ Pt plane. For that we shall 
calculate in the sequel the mean values of the action 5", the space-like plaquette Ps = Pxy 
and the transverse-like plaquette Pt = These quantities serve as order parameters 

which help us identify the various phases. 

A first attempt towards the determination of the phase diagram is through the mean 
field approximation. One considers the free energy in the mean field approximation, which 
(up to additive constants) is given by the expression: 

F{x) = -{(3s + 2(3t)[u'{x)]^ - u{x) + xu'{x). 

The function u{x) is defined through the relation 

exp[u(a;)] = + = 2cosh(x) u{x) = log(2) + log[cosh(x)] u\x) = tanh(a;). 

We observe that F{x) depends on the combination B = (3s + '2(3t- The free energy has 
always a local minimum at x = 0. For small B this is also the global minimum. As B 
increases, a second minimum shows up, which eventually wins and becomes the global 
minimum at B = 0.688. Thus the phase transition line is given by: 

^ 0.688 - (3s 
(3t= ^ . 

More accurately it is the segment of this line which corresponds to positive values of Ps 
and Pt- For Ps = we predict a critical value 0.344 for Pt, while as Ps increases the 
critical value decreases and finally becomes zero at Ps = 0.688. 

For the first set of measurements we have fixed Ps to several values, let Pt run and 
found the hysteresis loops which have been formed. The results of these measurements 
are displayed in figure 0. The subfigures correspond to the lattice volumes 6^, 8^, 10^ and 
12^ respectively. The line segments indicate the extents of the hysteresis loops. We have 
proceeded with steps of 0.005 for Pt and performed 200 iterations at each point. We have 
used plain Metropolis Monte Carlo as a simulation technique. The phase transition line 
tends towards the horizontal axis for large /?5. It is not clear from such measurements what 
happens for Ps —>■ oo, that is, whether the phase transition line meets the horizontal axis 
or it ends at some point. However, it is known analytically about the Pt = case that 
no phase transition should show up. Conceivably the transition weakens and eventually 
becomes a crossover before it meets the axis. 
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Figure 2: Phase diagram for (3s fixed at several values for 6'^, 8'^, 10'^, 12^ lattices. 

From the phase diagrams we may infer that the isotropic model will have a phase 
transition at = Pt ^ 0.25. This value is approximately the point where the Ps = Pt 
line meets the phase transition line. 

After this first overview of the phase structure we will study some of its characteristics 
in more detail. To this end we have performed long runs, sticking to a particular point 
of the parameter space, that is specific values for j3s and Pt and performing about three 
million iterations. When the parameters are near a first order phase transition we expect 
to see the eventual two-state signals. As a by-product of this procedure, the critical points 
may be determined with much greater accuracy than the one provided by the hysteresis 
loop method. We have found that, if we fix Ps, the critical Pt increases slightly with the 
volume, while for the isotropic model the volume dependence is very small. 

We now proceed with the presentation of the behaviour of the system in the immediate 
vicinity of the phase transition. The computer time evolution will be depicted for 8^ 
lattices. We concentrated on three values of Ps, namely 0.00, 0.20 and 0.60 and tried 
several values of Pt near the critical point. 

In figure ^ we show the time evolution of the transverse-like plaquette for a 8^ lattice 
for Ps = 0.00 and various values of Pt around 0.34. The system fiuctuates violently and, 
as Pt increases, it gradually spends its time more and more in the low Pt region. This 
behaviour provides evidence that the transition is of higher order. 

In figure ^jwe present the time evolution of the transverse plaquette for Ps = 0.20 in the 
phase transition region, that is for Pt around 0.26. We may clearly see the oscillation of the 
mean values between the two metastable states and the gradually increasing importance 
of the "low Pt" metastable phase with respect to the "high Pt" one as Pt increases. More 
precisely, the system starts by spending most of its "time" in the state with large Pt, but 
it gradually starts visiting also the state with the small Pt, until at some point it spends 
most of its time in the small Pt, as shown in the last subfigure. A very similar behaviour 
shows up for Ps = 0.30. 

For Ps = 0.60 the picture is strongly reminiscent of the Ps = 0.00 case. Figure ^ 
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shows relatively long runs for several values of I3t around 0.11. The system performs large 
fluctuations again and its mean value drifts toward small values of Pt as (3t increases. 

It appears that the phase transition is first order around the value (3s — 0.25, it 
remains strong for fis not very different from this value, but it weakens substantially for 
too small or too large values. 

Finally, in figure ^ we present the time evolution of the plaquette for the isotropic 
model on a 8^ lattice. For (3 around 0.24 we observe the phase transition and we may see 
the two metastable states. The fluctuations of the system between the two metastable 
states are very similar to the ones of the anisotropic model with (3 = 0.20. It appears that 
the phase transition is quite strong here, in agreement with previously obtained results. 

The final statement about the order of the phase transitions should come from a study 
of the volume dependence of the susceptibilities and the Binder cumulants. However, as 
one increases the volume, the system sticks to either of the metastable states and one 
cannot really observe the oscillation between the two states in a reasonable time. The 
only exception occurs for rather small lattice volumes. We have already presented the 
results for 8^ lattices, but it is difficult to observe something similar for larger volumes. 
However, one can easily see that if the susceptibility varies linearly with the volume 
(which is the sign of a first order transition), the gap between the two metastable states 
should be volume independent. Thus, we may get an idea about the order of the phase 
transitions by studying the volume dependence of the gap. If it is volume independent, 
we have a first order transition. If it decreases with the volume, there is a weaker phase 
transition (second or higher order). We find out that the gap does not actually depend 
on the volume. Figure |^ shows the results for Ps = 0.20. Thus it appears that this phase 
transition is of first order. The same picture also appears for Ps = 0.30. This means 
that there are first order phase transitions in a region around 0.25. For (3s away from this 
region one cannot really define a gap (compare figure ^). It is clear, however, that the 
phase transition is weak in this regime. 
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Figure 6: Long runs for the isotropic model at (3g = 0.236, 0.238, 0.240 and 0.242. 
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Figure 7: Tiie gap versus the volume. 
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